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Abstract 

We investigate systems of self-propelled particles with alignment interaction. 
Compared to previous work [13^ I18j . the force acting on the particles is not nor- 
malized and this modification gives rise to phase transitions from disordered states 
at low density to aligned states at high densities. This model is the space inho- 
mogeneous extension of [19] in which the existence and stability of the equilibrium 
states were investigated. When the density is lower than a threshold value, the 
dynamics is described by a non-linear diffusion equation. By contrast, when the 
density is larger than this threshold value, the dynamics is described by a similar 
hydrodynamic model for self-alignment interactions as derived in |131ll8j. However, 
the modified normalization of the force gives rise to different convection speeds and 
the resulting model may lose its hyperbolicity in some regions of the state space. 
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1 Introduction 



The context of this paper is that of [I3] and is concerned with a kinetic model for self- 
propelled particles and its hydrodynamic or diffusion limits. The particles move with 
the same constant speed and their velocity directions (which belong to the sphere S) 
align to the local average orientation, up to the addition of some noise. This model has 
been proposed as a variant of the Vicsek particle model [35]. In this paper, we remove 
the normalization of the force intensity which was done in [13]. This apparently minor 
modification leads to the appearance of phase transitions, which have been studied in 
the space- homogeneous setting in [TH]. In [TH], it is proved that the equilibrium distribu- 
tion function changes type according to whether the density is below or above a certain 
threshold value. Below this value, the only equilibrium distribution is isotropic in ve- 
locity direction and is stable. Any initial distribution relaxes exponentially fast to this 
isotropic equilibrium state. By contrast, when the density is above the threshold, a sec- 
ond class of anisotropic equilibria formed by Von-Mises-Fischer distributions of arbitrary 
orientation appears. The isotropic equilibria become unstable and any initial distribution 
relaxes towards one of these anisotropic states with exponential speed of convergence. 
We would like to emphasize the connection of the presented alignment models to the the 
Doi-Onsager [T6| [27] and Maier-Saupe [23] models for phase transition in polymers. The 
occurrence of phase transitions makes a strong difference in the resulting macroscopic 
models as compared with the ones found in [131 HH] , where no such phase transitions were 
present. 

In the present paper, we rely on this previous analysis to study the large-scale limit of 
the space-inhomogeneous system. In the regions where the density is below the threshold, 
the convection speed becomes zero and the large-scale dynamics becomes a nonlinear 
diffusion. On the other hand, in the region where the density is above the threshold, 
the large-scale dynamics is described by a similar hydrodynamic model for self- alignment 
interactions as derived in [T3| [T8]. However, the modified normalization of the force gives 
rise to different convection speeds and the resulting model may lose its hyperbolicity in 
some regions of the state space. 

The Vicsek model [HS], among other phenomena, models the behaviour of individuals 
in animal groups such as fish schools, bird flocks, herds of mammalians, etc (see also 
[H El El EO]). This particle model (also called 'Individual-Based Model' or 'Agent-Based 
model') consists of a discrete stochastic system for the particle positions and velocities. 
A time-continuous version of the Vicsek model and its kinetic formulation have been 
proposed in [13]. The rigorous derivation of this kinetic model has been performed in [1]. 

Hydrodynamic models are more efficient than particle models for large numbers of 
particles, because they simply encode the different particles quantities into simple av- 
erages, such as the density or mean-velocity. We refer to [3 [T71 [Ml [251 1321 E31 [31] for 
other models of self-propelled particle interactions. Rigorous derivations of hydrodynamic 
models from kinetic ones for self-propelled particles are scarce and [T3l [T8] are among the 
first ones (see also some phenomenological derivations in fI2\ [2H1 [30])- Similar models 
have also been found in relation to the so-called Persistent Turning Walker model of fish 
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behavior [T21II1]. Diffusive corrections have also been computed in [15]. We refer to [SlEj 
for other macroscopic models of swarming particle systems derived from kinetic theory. 
In particular, we mention [18] where a vision angle and the dependence of alignment 
frequency upon local density have been investigated. 

The outline of this paper is as follows. In Section |2l we describe the Individual- 
Based Model (IBM), and its mean-field limit. In Section [SI we investigate the properties 
of the rescaled mean-field model. We prove that there are two possibilities for a local 
equilibrium, depending on the value of its density p. Section H] is devoted to the derivation 
of the diffusion model when the density p is below the threshold. Finally, in Section [5l 
we derive the hydrodynamic model for self-alignment interactions in the region where p is 
above the threshold and study its hyperbolicity. A conclusion is drawn in section [HI Two 
appendices are added. In appendix 1, we calculate a Poincare constant which provides us 
with a fine estimate of the rate of convergence to the equilibrium states. In appendix 2. 
some numerical computations of the coefficients of the model are given. 



2 Particle system and mean-field limit 

We consider oriented particles in M", described by their positions Xi, . . . Xi^ and their 
orientation vectors wi, . . . , ujn belonging to S, the unit sphere of M". We define the mean 
momentum of the neighbors of the particle k by 

1 ^ 
i=i 

In this paper, the observation kernel K will be supposed isotropic (depending only on 
the distance \Xj — X^l between the particle and its neighbors), smooth and with compact 
support. Introducing a non-isotropic observation kernel, as in [TSj would lead to the same 
conclusion, with a slightly different convection speed for the orientation in the macroscopic 
model, but the computations are more complicated. Therefore we focus on an isotropic 
observation kernel for the sake of simplicity. 

The particles satisfy the following system of coupled stochastic differential equations 
(which must be understood in the Stratonovich sense), for k G |1,A^]: 

dXk = uJk^t (2.1) 
duk = (Id - cuk ® cok)Jk dt + V2d{ld -ujk® uJk) o d5f , (2.2) 

The first equation expresses the fact that particles move at constant speed equal to unity, 
following their orientation Uk- The terms stand for N independent standard Brownian 
motions on M", and the projection term (Id — tUfc ® Uk) (projection orthogonally to Uk) 
constrains the norm of Uk to be 1. We have that (Id — cufc ® 0Jk)Jk = ^uj{(^ ■ </fc)U=ajfe, 
where V^^ is the tangential gradient on the sphere. So the second equation can be under- 
stood as a relaxation (with a rate proportional to the norm of Jk) towards a unit vector 
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in the direction of J^, subjected to a Brownian motion on the sphere with intensity -\/2rf- 
We refer to for more details on Brownian motions on Riemannian manifolds. 

The interaction term (first term of (12.21) ) is the sum of smooth binary interactions. 
This model is an intermediate between the Cucker-Smale model [9], where there is no 
constraint on the velocity and no noise, and the time-continuous version of the Vicsek 
model proposed in [13], where the velocity is constant and noise is added. Indeed, in 
[T^ . Jfc is replaced by z/f2fc, where ilk = is the unit vector in the direction of Jk and 
the relaxation frequency z/ is a constant. Therefore, in [13], the interaction term cannot 
be recast as a sum of binary interactions and has a singularity when Jk is close to 0. 
The model presented here brings a modification consisting in letting u depend (linearly) 
on the norm of the velocity Jk- A related modification has previously been introduced 
in [IH], consisting in letting the relaxation parameter u depend on a local density pk, but 
the modification considered here brings newer phase transition phenomena. 

From the Individual-Based Model (12. ip . (12. 2p . we derive a mean- field limit as the 
number of particles tends to infinity. We define the empirical distribution by 

1 ^ 

f^{x,uj,t) = — 5]]5(x,w,a;,(t))(a;,w), 
1=1 

where the Dirac distribution is defined by duality by v)k"xs = '^{^,^) for any 

smooth function (p G C(W^ x S), the duality product (■, ■)Rnxs extending the usual inner 
product of L^(]R" x S). For convenience, the integration measure is supposed of total mass 
equal to 1 on the sphere S, and we have (/^, l)Mnxs = 1- Denoting the convolution with 
respect to the space variable by *, and the duality product on the sphere by (-, ■)§, we 
get Jk = {K * f'^{Xk),co)n- If there is no noise {d = 0), it is easy to see that satisfies 
the following partial differential equation (in the sense of distributions): 

dtf'' + u ■ V./^ + V. ■ ((Id -u(g) uj)JfNf) = 0, 

where V^j- denotes the divergence operator on the unit sphere, and 

JfN{x,t) = {{K*f){x),uj)s. 

When noise is present (rf 7^ 0), the empirical distribution tends to a probability density 
function / satisfying the following partial differential equation: 

dtf + uj- VJ + ■ ((Id -u^ u)Jff) = dAJ, (2.3) 

with 

Jf{x,t)= {K*f){x,u,t)udu. (2.4) 
Js 

This result has been shown in [1], under the assumption that the kernel K is Lipschitz 
and bounded. 

Eqs. (12. 3p . (12. 4p are the starting point of our study. We notice that there is a com- 
petition between the alignment and diffusion terms. The alignment term is quadratic 
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while the diffusion term is hnear. So we can expect that ahgnment wins over diffusion 
for high densities while at low densities, diffusion dominates. This is the source of the 
phase transition rigorously studied in the space-homogeneous setting in [19]. In this ref- 
erence, it is proven that there is a unique isotropic equilibrium at low density but beyond 
a certain density threshold, another family of non-isotropic equilibria in the form of Von- 
Mises-Fischer distributions arises. Above this threshold, the isotropic equilibria become 
unstable and the anisotropic ones become the stable ones. Therefore, we expect different 
large-scale limits according to whether the density is lower or larger than this threshold. 

We now make some preliminary remarks and assumptions. We suppose that the 
kernel K is integrable, and that its total weight Kq = K{x)dx is positive. Writing 

f{x,u,t) = f{^x,uj,^t) and K{x) = j^K{^x), 
we get that / satisfies ( 12. 3 p with d = 1 and K replaced by K in (12. 4p . and we have 




K{x)dx = 1. 



So without loss of generality, we can suppose that d = 1 and that Kq = 1. 

We are now ready to investigate the large-scale behavior of (12. 3p . ( 12.41) in space and 
time. The derivation of the macroscopic limit proceeds as in [13], and follows closely the 
presentation of [H] , so we only give a summary, focusing on the points which are specific 
to the present model, in particular the distinction between the ordered and disordered 
phases. 



3 The macroscopic limit 
3.1 Hydrodynamic scaling 

In order to observe the system at large scales, we perform a hydrodynamic scaling. 
We introduce a small parameter e, and the change of variables x' = ex, t' = et. We 
write f'^{x',uj,t') = f{x,u,t), and K^{x') = -^Ki^x). Then satisfies 

e{dtf + u ■ V.r) = -V. ■ ((Id - a; ® a;)^^) + A^f , (3.1) 

with 

J^x,t)= l{K'*f){x,u,t)uduj. (3.2) 
Js 

The purpose of this paper is to derive a formal limit of this rescaled mean-field model 
when the parameter e tends to 0. The first effect of this hydrodynamic scaling is that, up 
to order 1 in e, the equation becomes local. Indeed, supposing that does not present 
any pathological behavior as e — 0, we get the following expansion: 

J},{t,x) = Jf.{t,x) + 0{e^), (3.3) 
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where the local flux Jf is defined by 



Jf{x,t) = / f{x,u,t)udu. (3.4) 

The proof of this expansion is elementary and omitted here (see e.g. Appendix A.l of [18j). 
We also define the density pf associated to / by 



pf{x,t) = / f{x,uj,t)duj. (3.5) 
J§ 

Hence, Eq. 03. ip becomes, after dropping the 0{e'^) term: 

s{dtr + u-VJ') = Q{r), (3.6) 

with 

Qif) = -V. ■ ((Id -uj(S) cu)Jff) + A^f. (3.7) 

This paper is concerned with the formal limit e — )■ of this problem. 

We remark that the collision operator Q acts on the u variable only. The derivation 
of the macroscopic model relies on the properties of this operator. An obvious remark is 
that 



[ Q{f)doj = Q (3.^ 



which expresses the local conservation of mass. The first step of the study consists 
in characterizing the equilibria, i.e. the functions / such that Q{f) = 0. Indeed, 
when e — )■ 0, Qif^) — )■ and the limit / = limg^o belongs to the set of equilibria. 
This characterization is the purpose of the next subsection. 

3.2 Equilibria 

For any unit vector f2 G §, and k ^ 0, we define the so-called Von-Mises-Fisher distribu- 
tion [SH] with concentration parameter k, and orientation Q by 

We note that the denominator depends only on k. Mi^q is a probability density on the 
sphere, and we will denote by (■)A/Kn the average over this probability measure. For 
functions 7 depending only on u -Q, the average (7(0; ■ Q))m^q does not depend on fl and 
will be denoted by (7(cos 6'))^/^. Using spherical coordinates, this average is given by: 

£ 7(cos 6) e^=°^ ^ sin"-2 6 dO 
The fiux of the Von-Mises-Fisher distribution is given by 

JM^n = {^)M,n = c{k)Q, (3.10) 
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where the so-called order parameter c(k) is such that ^ c(k) ^ 1 and is defined by 

/o"cos^e'''=°^''sin'^-2ed^ 



c{k) = {cos6)m^ 



(3.11) 



c(k) measures how the distribution M^^ is concentrated about Q. When c(k) = 0, M^n 
is the uniform distribution M^t^ = 1, and when c{k) — j- 1, we have M^t^ — > S^{uj). 

We remark that the dependence of M^n upon k, and Q only appears through the 
product kQ. In this way, we can consider Mj for any given vector J G M". We also note 
that Va;(Mj) = (Id - a; (g) uj)J Mj. Therefore 



QU) = V. ■ 
Using Green's formula, we have 



Is 



Qif) 
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Mj, 



duo 



f 



f 



Mj, 



M, 



■ V„ 



M, 



Mj, du. 



and 



Is 



Qif) 



f 



M 



duo 



f 



M 



Mj, duo ^ 0. 



(3.12) 



Definition 3.1. A function f{u}) is said to be an equilibrium of Q if and only if Q{f) = 0. 



Let / be an equilibrium. Using fl3.12p . we deduce that jj— is a constant. There- 

fore, / = pfMjj is of the form pMi^Q with k ^ and 1) G S (we note that in the 
case \Jf \ =0, then k = and we can take any f2 G S because / is then just the uniform 
distribution). Using fl3.10p . we get 

nQ = Jf = pJm^q — pc{f^)^, 

which leads to the following equation for k (compatibility condition): 

Pc(k) = K. (3.13) 

The study of this condition and the classification of the equilibria can be found in [TU] . 
The key point is to notice that the function k is decreasing and tends to ^ as k — )■ 0. 

Therefore, there is no other solution than k, = if p ^ n. By contrast, if p > n, there is 
a unique strictly positive solution in addition to the trivial solution k, = 0. This leads to 
the following proposition. 



7 



Proposition 3.2. (i) If p ^ n, n = is the only solution to the compatibility relation 
\3.1'J\) . The only equilibria are the isotropic ones f = p, with arbitrary p ^ 0. 

(a) If p > n, the compatibility relation ^3.13^ has exactly two roots: k = and a unique 
strictly positive root denoted by k{p). The set of equilibria associated to the root k = 
consists of the isotropic equilibria f = p, with arbitrary p > n. The set of equilibria 
associated to the root /t(p) consist of the Von Mises-Fischer distributions pM^i^p^Q^ with 
arbitrary p > n and arbitrary f2 G § and forms a manifold of dimension n. 

The rate of convergence to the equihbria have been studied in [19] in the spatially 
homogeneous setting. The results are recalled in the next section. 

3.3 Rates of convergence to equilibrium in the spatially homo- 
geneous setting 

Denoting by g"^ = f^/pf^ the velocity probability distribution function, we can rewrite 
(13.61) under the following form (omitting the superscripts e for the sake of clarity and 
neglecting the 0{e'^) term): 

e{dt{pg) + u ■ V,{pg)) = -(p)'V^ ■ ((Id - w ® u)Jgg) + pA^g. 

In the spatially homogeneous setting, we let Vx{pg) = and get 

edtipg) = -{pfV^ ■ ((Id -uj® u)Jgg) + pA^g = Q{pg). (3.14) 

Integrating this equation with respect to u and using (13. 8p . we find that dtp = 0. There- 
fore, p is independent of t and can be cancelled out. The homogeneous equation (I3.14p 
therefore takes the form: 

edtg = -pV^- ((Id - oj ® uj)Jgg) + A^g. (3.15) 

We now remind the definitions of global and asymptotic rate. 

Definition 3.3. Let X be a Banach space with norm || ■ || and let f{t): R+ X be a 
function oft with values in X . We say that f{t) converges to foo with global rate r if and 
only if there exists a constant C which only depends on ||/o|| , such that 

||/(t)-/ooKCe-^*. (3.16) 

We say that f[t) converges to f^o with asymptotic rate r if and only if there exists a 
constant C depending on /o (but not only on \\fo\\) such that Ii3.16]) holds. Finally, we 
say that f(t) converges to foo with asymptotic algebraic rate a if and only if there exists 
a constant C depending on fo 

\\f{t)-U\^^. 

Now, concerning problem (13.151) . we can state the following theorem: 
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Theorem 3.4. |19] Suppose Qq is a probability measure, belonging to H^{S). There exists 
a unique weak solution g to /i3.15\) . with initial condition g{0) = g^. Furthermore, this 
solution is a classical one, is positive for all time t > 0, and belongs to C°°((0, +00) x S). 

(i) If Jg^^ 7^ 0, the large time behavior of the solution is given by one of the three cases 
below: 

- Case p <n: g converges exponentially fast to the uniform distribution, with global rate 

ne 

in any norm. 

- Case p > n: There exists f2 G S such that g converges exponentially fast to Mk(p)q, with 
asymptotic rate greater than 

r{p) = ^[pc(k(p))2 + n- p]A,(p) > 0, 
in any norm, where is the best constant for the following Poincare inequality: 

(IV^Hm.. > - {9)M.nf)M.a, (3-18) 

We have 

r(p) ~ -2(n - - 1), when p ^ n. (3.19) 
e n 

- Case p = n: then g converges to the uniform distribution in any norm, with algebraic 
asymptotic rate 1/2. More precisely, we have: 



\g-nHv<C ( 



t. 



(a) If J go = 0: Then, 1^3. 15\} reduces to the heat equation on the sphere. So g converges 
to the uniform distribution, exponentially fast, with global rate r = ^ in any norm. 

Remark 3.1. That go is a probability measure implies that go G H'^{E>) for all s < 
However, the theorem holds for all s. So for s ^ that go G H^(E>) is not a mere 

consequence of being a probability measure and must added to the hypothesis. 

Now, we comment the resuhs of this theorem. First, in the supercritical case (when p > 
n), the uniform distribution is an unstable equilibrium: for any perturbation g of the 
uniform distribution such that Jg ^ 0, the associated solution converges to a given Von- 
Mises distribution, with a fixed concentration parameter k(p) defined by the compatibility 
condition fl3.13p . Second, the rates of convergence to the equilibrium are exponential. In 
the supercritical case, these rates are only asymptotic ones, but we can prove a uniform 
bound on these rates for p in any compact interval. A more precise study of the behavior 
of these rates is left to future work. 
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Figure 1: Rates of convergence to equilibria in dimensions 2, 3, and 4, as functions of the 
density p. 

Therefore, when e is small, the function converges rapidly to a given equilibrium, 
provided that the rate satisfies r(p) — )■ oo when e — )■ 0. In the case p < n, thanks to 
f l3.17p . this condition is equivalent to saying that e = o{n — p). In the case p > n, thanks 
to f lS.igp . the condition e = o{n — p) implies that r(p) — )■ oo when e — )■ uniformly in any 
bounded p interval of the form [n, A] with A < oo. However, a uniform estimate from 
below of r{p) is lacking when p — )■ oo. But we can reasonably conjecture that away from 
a buffer region \p — n\ = 0{e), the convergence to the equilibrium is exponentially fast. 

Some elements towards a uniform estimate of the rate r{p) are provided in Appendix 
1. Furthermore, in Appendix 2, we compute A^ and then r{p) numerically. The results 
are depicted in Fig. [T]for dimensions 2, 3, and 4. We observe that for p > n, r{p) grows 
linearly with p, which supports our conjecture. 

Therefore, in the general space-inhomogeneous case, we will assume that the formal 
limit of as e — )■ is given by a function f{x,u},t) which has a different velocity profile 
according to the position of the local density p{x, t) with respect to the threshold value n. 
For this purpose, we define the disordered region TZd and the ordered region TZo as 



T^d = {{x^t) I ^ P^{^^ t) ^ e, as £ — 7- }, 
TZo = {(x, t) I p^(x, t) — e, as e — )• }. 



(3.20) 
(3.21) 



We assume that as e — )■ we have 



/"(x, w, t) p{x, t), V(x, t) e 7^d, 

f{x, u, t) p(x, t) M«(p)Q(^,t), V(x, t) E n, 




(3.23) 



and that the convergence is as smooth as needed. 
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The goal is now to derive evolution equations for p{x,t) and Q{x,t). This is the 
subject of the following two sections. We already note that, integrating fl3.6p with respect 
to u and using (13. Sp . we get the mass conservation equation 

dtp' + V,-iJf.) = 0. (3.24) 

4 Diffusion model in the disordered region 

We derive the macroscopic model in the disordered region TZd C M"^, using (I3.22p . With (13.240 
and the fact that J/c — )• = 0, the mass conservation equation reduces to 

dtp = 0. 

To obtain more precise information, we look for the next order in e, using a Chapman- 
Enskog method, similarly to the case of rarefied gas dynamics (see [10] for a review). We 
prove the following theorem: 

Theorem 4.1. When e tends to zero, the (formal) first order approximation to the solu- 
tion of the rescaled mean-field system liS. 6\) . IjjS. 7| j in the disordered region TZd defined by 
I13.2(J\) is given by 

^«-/ X «■/ X nu ■ VxP'^ix.t) , 

r{x, u, t) = X, t) - e \ \ 4.1 

[n — l)(n — p^(x, t)) 

where the density p' satisfies the following diffusion equation 

dtp' = ^ V, ■ Vxp') . (4.2) 

n — 1 \n — p^ J 



Proof. We let p' = pfe and write f = p'{x,t) + efl{x,u,t) with J^ffdu = 0. Inserting 
this Ansatz into (13. 4p . we get 

and the model (13. 6p . (13. 7p becomes: 

dtp' + u ■ Va^p' + e{dt + u ■ Vx)fl = -V^((Id - w ® w) J/f p") + A^/f 

~ eV^{{\d-uj®uj)Jfep'). 



Additionally, gives: 

dtp' + eW.,-{Jp^) = Q. (4.4) 

In particular dtp' = 0{e). We need to compute to find the expression of the current. 
But, with this aim, we may retain only the terms of order in (14.30 . Since 

Va;((Id -oj® uj)A) = -{n - 1)A ■ oo, 



11 



for any constant vector A G M", the equation for ff reads: 

A^f! = (V,p^ -in- 1)/ J^j) ■ u + 0{e). 

This equation can be easily solved, since the right-hand side is a spherical harmonic of 
degree 1 (i.e. is of the form A-u; we recall that A^^^A -u) = —{n — 1)A ■ u and that A ■ u 
is of zero mean). Then: 

f! = -;^(V./ -in- l)p^Jf.) ■ uj + 0{e). 

We immediately deduce, using that w ® ujdu = ^Id: 

which implies that 

Inserting this equation into (14. 4p leads to the diffusion model (14.21) and ends the proof. ■ 

Remark 4.1. The expression of ff, which is given by the 0{e) term of ( [^. j| ) confirms that 
the approximation is only valid when n — p"^ ^ e. The diffusion coefficient is only positive 
in the disordered region and it blows up as tends to n, showing that the Chapman- Enskog 
expansion loses its validity. 

5 Hydro dynamic model in the ordered region 
5.1 Derivation of the model 

We now turn to the ordered region TZo C defined by (13.211) . The purpose of this section 
is to give a formal proof of the following: 

Theorem 5.1. When e tends to zero, the (formal) limit to the solution /'^(x,Co',t) of the 
rescaled mean-field system liS. 6\) . IjjS. 7] ), in the ordered region TZq C M" defined by Ii3.21\) . 



is given by 

fix,u,t) = p{x,t) M^^p^.,^t))n{x,t){(^), (5.1) 

where the Von- Mis es- Fischer distribution M^q_ is defined at (13. 9p . and the parameter n 
is the unique positive solution to the compatibility condition (I3.13p . Moreover, the den- 
sity p > n and the orientation Q E S satisfy the following system of first order partial 
differential equations: 

dtp + V, ■ (pen) = 0, (5.2) 

p{dtn + c{n ■ v^)n) + A(id -n^ i^) v^p = o, (5.3) 
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where the coefficient c = c{k{p)) is defined at fl3.11l) . the coefficient c = c{k{p)) will he 
defined later on at fl5.9p . and the parameter A = A(p) is given by 

. p-n- Kc ^, 

A = -f -. 5.4 

K[p — n — KC) 

Proof: From now on, we will drop the dependence on p in the coefficients when no 
confusion is possible. With f l3.23p . ^ f, where / is the stable local equilibrium (15. ip . 
We now derive the evolution equations (15. 2p . (15. 3 p for p and Q. 

We recall that the concentration parameter k, satisfies the compatibility condition (I3.13P 
where the order parameter c is defined by (13.111) and that we have Jj = pcVt. Therefore, 
eq. (I3.24P in the limit e — )■ 0, reads 

dtp + V, ■ {pen) = 0. 

To compute the evolution equation for Q, the method proposed originally in [13] 
consists in introducing the notion of generalized coUisional invariant (GCI). This method 
has been then applied to [HI HE]- The first step is the definition and determination of the 
GCI's. We define the linear operator L^n associated to a concentration parameter k and 
a direction Q as follows: 



L^nif) = A^/ - kV^ ■ ((Id -cu® Lo)nf) = V, 



/ 



M, 



so that Q{f) = Lj^.{f). We define the set C^n of GCI's associated to k G M and Q E S 
by: 

C.n = (V'l / L^n{f)4^duj = 0, V/ such that {Id - Q ^ Q)Jf = 1 . 
Hence, if ijj is a GCI associated to k and Q, we have: 

/ Qif) V^dcu = 0, V/ such that J/ = kH. 



The determination of C^n closely follows [18]. We define the space 

V = {g\{n- 2)(sin^)t"2^ e L\0, vr), (sin^)t"i(^ ^ ^1(q^ ^^y^ ^^^^ 

and we denote by g,^ the unique solution in V of the elliptic problem 

L:gie) = sine, (5.6) 

where 

Llgie) = -{sinef-e—'^{isiner-'e^-'^{e)) + i^,g{e). (5.7) 
Then defining /i^ by g^^O) = hf^{cos6) sin 6', we get 

C^n = {h^ioo ■n)A-oj + C \ C G M, AeW, with A-n = 0}. 
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Ci^Q is a vector space of dimension n, since A is a vector with n—1 independent components. 

The next step consists in multiplying f l3.6p by a GCI associated to and fi^ such 
that Jfs = k'^Q'^, and to integrate it with respect to u. For any vector A G M", with A-Q"^ = 
0, we get 

/ QiDK^iuj -n') A-codco = 0. 



So, the vector 



is parallel to or equivalently (Id — fi^ ® = 0. Using (13.61) . we get: 

In the limit e — )■ 0, we get 

{ld-n(»n)X = 0, (5. 

where 



X= / (9t(pM,n)+c^- V,(pM,n))/i«(wf^)u;da;. 
Finally it has been proved in [TH] that (15. 8 p is equivalent to (15. 3 p with 



/J'cos^/i«(cos^)e^'=°^^ sin"^d^ 
/Q"/i^(cos^)e^^°«^ sin"^d^ 



c = {cos 9)^^ = ^° , (5.9) 



1 p d/t , , , 

A = - + ^— ^-c. 5.10 
K n dp 

We can now compute a simpler expression of A. We differentiate the compatibility 
condition (I3.13P with respect to k, and we get 

dp dc ^ 
dn ^dn 



We have 



dc _ d //o"cos^e^=°'^''sin"^2^d^ 



dn dK\ /Q"e«^°«^sin"-2ed^ 

cos2 ^ e''™^'' sin"~2 0^9 f cos ^ e'^™'^^ sin""^ ^ 



/J" e'^ '^"^ « sin'^-^ ^ d^ V lo e'^ ^ sin"-^ 9 d9 
/J'sin2^e^^°^^sin""2^d^ 2 
re'^^°«'^sin"-2^d^ ^ 



c 



1 - (n - 1)- - c- 



2 
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Therefore we get 

dp K, dp dc / / N c 9, , , 

= - ^ = I - p— = I - p{l - {n - I)- - c^) = n - p + Kc, 5.11 
QK p dK (Ik k. 

and finally 

1 c — c n — p + Kc 

^ = - + = —r ^' 

K n — p + KC K[n — p + KC) 

which ends the proof of Theorem 15. 1[ ■ 

The next part is devoted to the study of the properties of the model f l5.2p -( l5r3|) in the 
ordered region. 



5.2 Hyperbolicity of the hydrodynamic model in the ordered 
region 

We first investigate the hyperbolicity of the hydrodynamic model (15. 2p - ( 15731) . We recall 
some definitions. Let 

n 

dtU + J2MU)d..U = 0, (5.12) 

i=l 

be a first order system where x G M", t ^ 0, U = (f/i, . . . Um) is a m-dimensional vector 
and {Ai{U))i=i^,,,^n are n m x m-dimensional matrices. Let Uo G M™. The constant and 
uniform state U{x,t) = ?7o is a particular solution of (I5.12p . The linearization of (I5.12p 
about this constant and uniform state leads to the following linearized system: 

n 

dtU + J2MUo)d^.u = 0. (5.13) 

i=l 

We look for solutions of (I5.13P in the form of plane waves u{x, t) = u e*('^'^'"'^*), with A; G 
and a; G C. Such solutions exist if and only if uj/\k\ is an eigenvalue of the matrix A{k/\k\) 
and u is the related eigenvector, where for a direction ^ G S, the matrix A(^) is defined 
by 

n 

A{0 = J2A{U)^,. (5.14) 

1=1 

The problem (I5.12p is said to be hyperbolic about Uo, if only purely propagative plane 
waves with real u can exist or equivalently, if A{C,) has real eigenvalues for any ^. We also 
must rule out polynomially increasing in time solutions which could exist if the matrix 
would not be diagonalizable. This leads to the following definitions: 

Definition 5.2. (i) Let Uq G M™. System Ii5.12\) is hyperbolic about Uq if and only if for 
all directions G S, the matrix A{C,) is diagonalizable with real eigenvalues. 

(a) System 1^5. is hyperbolic, if and only if it is hyperbolic about any state Uq in the 
domain of definition of the matrices Ai{U). 
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The linearization of system f l5.2p - fl5.3p about a stationary uniform state {po,Qo) is 
obtained by inserting the following expansion 



p = Po + 6r + o{6), (5.15) 
n = no + 6W + o{d), (5.16) 

with 5 <^ 1 a small parameter and r = r(x, t), W = W{x, t), the first order perturbations 
of p and Q. Given that \Q\ = \Qq\ = 1, we have W ■ Qq = 0. Inserting 05.151) . fl5.16p 
into fl5.2l) -( !531) leads to the following linearized system: 



dtr + 7o(l]o ■ V.)r + poCo(V. ■W) = 0, (5.17) 

dtW + co(fio ■ V.)Vr + —(Id - 1^0 ® f^o) V.r = 0, (5.18) 

Po 

W-no = 0, (5.19) 



with 

dc 

7(p) = c + p^, 

and 7o = 7(po), Cq = c(po), Cq = c(po ) and Ap = A(po)- 

Next, we show that system f l5.17p -( !5.19p is invariant under rotations. This will allow us 
to choose one arbitrary direction C, in the definition f l5.14p instead of checking all possible 
directions. For this purpose, let i? be a rotation matrix of M", i.e. Ris anxn matrix such 
that = R~^, where the exponent T denotes transposition. We introduce the change 
of variables x = Rx' and define new unknowns 

r{x) = r'(x'), W{x) = RW'{x'), = R^'^. 

We note the following identities 

n'^ ■ W'{x') = no ■ W{x) = 0, 
V,r(x) = RV,>r'{x'), 
V^W{x) = RV^W{x')R^, 
{V^-W){x) = {V,>-W'){x'), 

(fio ■ V,.)W^(x) = iV^Wix)fno = R{V^Wix')f% = Rin'o- V^>)W'{x'), 
(fio- V.)r(x) = ((][,■ V.Or'(x'). 

With these identities, it is easy to show that (r', W) satisfies system f l5.17p -( !5rT9l) with Qq 
replaced by Q'q. 

The rotational invariance of fl5.17p - fl5.19p shows that, in order to check the hyperbolic- 
ity, it is enough to choose any particular direction ^. Let us call this arbitrary direction z, 
with unit vector in this direction denoted by e^. To check the hyperbolicity of waves 
propagating in the z direction it is sufficient to look at the system where all unknowns 
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Figure 2: Coefficient A in dimensions 2, 3, and 4. 



only depend only on the space coordinate z and on the time t. Denoting by 6 the angle 
between the z direction and Q, we can write: 

Q = cosOcz + sm9v, 6'G[0,7r], f G §„_2 5 

where S„,_2 is the sphere of dimension n — 2 collecting all unit vectors orthogonal to Cz- 
With these hypotheses, system fl5.2lj -( !531) is written. 



dtp+ dz{pc{p) cos 6 
p[dt{cos6) + c{p) cos 6 dz{cos 9)] 
dtv + c{p) cosOdzV = 0, with |t>| 



= 0. 

A sin^eOzp = 0. 
: 1 and ■ V = 0. 



(5.20) 
(5.21) 
(5.22) 



In the special case of dimension n = 2, the system reduces to fl5.20p - fl5.21l) . with 6 G 
(— 7r,7r) and Q = cosOe^ + sin6'wo, where vq is one of the two unit vectors orthogonal 
to Cz- 

The hyperbolicity of this system depends on the sign of A. Proposition 15.51 below 
shows that A < in the two limits p — n and p — ?■ oo. Additionally, the numerical 
computation of A, displayed in Fig. [21 provides evidence that A < for all values of p, at 
least in dimensions n = 2, 3, and 4. Therefore, we assume that 



A < 0. 

We first check the local hyperbolicity criterion: 



(5.23) 



Proposition 5.3. We assume / (5. 23\) . Then, system l{5.2IJ\) - [5. 2S\) is hyperbolic about (p, 6, v) 
if and only if 

|c ■ 



tan^l < tanOr 



n—p+KC I 



2^/^ 



(5.24) 
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Proof: We apply [18] and find that the hyperbohcity criterion is written: 

|c- j-{pc)\ 
Itan^l < !^==^. 

Using the compatibihty condition f l3.13p and (15. lip , eq. f l5.24p follows. ■ 

As for global hyperbohcity, we have 
Proposition 5.4. We assume /i5.2!^) . Then, system li5.2(J\) - l[5.2^) is not hyperbolic. 

Proof: It has been proved in [18] that system f l5.2p -( l573|) is hyperbolic if and only if A > 0. 
As we assume f l5.23p . it follows that the system is not hyperbolic. ■ 

We now provide asymptotic expansions of the coefficients which show that, at least 
when p — )■ n or p — )■ oo, we have A < 0. 

Proposition 5.5. We have the following expansions: 

(i) When p — )■ n: 

^=^v^ + 0(p-n). 



(^c = 1 - -Tri^Vp^ + 0{p - n) 



(a) When p — 7- oo.- 



C = 1 - ^p-' + (nzlM!l±l)p-2 + 0(p-3), 
C=l- «±lp-l _ ill±lM^p^2 _^ 0(^-3) 

X = -^p~' + Oip-% 



= arctan( ^^ ) + 0(p-^). 



Proof: Using the compatibility condition fl3.13p . the expression fl5.4p depends only on k, c, 
and c. With the asymptotic expansion of c and c as k — and k — )■ oo given in [TH] , we 
can get an expansion for A. We have 



c 



c 



- + as «: ^ 0, 

-^ + ^""8ir'^ +Q(^"') asK^oo, 



l-i^ + ^^l±^ + 0{.~^) as.^oo. 
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We first compute an expansion of p = -. We get 



^~^K + ^+^-^^^ + 0{^-') as«:^oo. ^^^^^ 



Using tlie definition (15 .4^ . we then get 



A 



-^ + 0(1) asK^O 
■|i^ + 0(K"^) as K oo. 



We can also expand the threshold angle 6^ in terms of n. We get 

Or 



\arctan( ^^ ) + 0{k~^) as k ^ oo. 



We can now reverse the expansion fl5.25p to get an expansion of k, (and then of the other 
coefficients) in terms of the density p. We get 



y/n + 2y/p — n + 0{p — n) as p — )■ n, 

t-l _ (rt-l)(n 
2 8p 



p-^-(^^^^ + 0(p-2) asp^oo. 



Inserting this expansion into the previous ones, we finally deduce the expressions stated 
in proposition 15.51 ■ 

When p ~ n, since |A| = —A is large compared to p, which is large compared to pc, 
the behavior of the orientation equation (I5.3p can be compared to the behavior of 

d^VL = ^ (Id - 1] ® VL)V.^p, 
P 

which relaxes Q to the unit vector VxP/\'VxP\, with rate 



P 



This actually makes sense only if the rate of convergence to the equilibrium ^r(p) ~ 
^^^(p — n) in the neighborhood of n is large compared to this relaxation rate. This 
requires 5 ^ (p — n)2\^^p\. In this case the leading behavior of the system is given by 

which is an ill-posed problem, being some kind of nonlinear backwards heat equation. 
To stabilize this system, a possibility is to derive a first order diffusive correction to 
model (15. 2p . (15. 3 p using a Chapman- Enskog expansion. Such a correction has been derived 
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Figure 3: The velocities c and c in dimension 2 (left) and 3 (right). 

in [15] for the model of [I3] , but leads to complicated terms. Another possibility is to add 
some contribution of the non-locality of the interaction in the spirit of [TTj . 

When p is large, c and c are close to 1, and A is small. In the intermediate regime, 
numerical computations (see Appendix 2) show that there is a significant difference be- 
tween c and c. This means that the information about velocity orientation travels slower 
than the fluid. Fig. [3] displays c and c as functions of p, in dimension 2 and 3. 

Finally, when p — )■ oo, the critical angle 9c tends to a positive value arctan( ^'^"^^^ ). 
Numerically, we see that 9c is always larger than this limit value. Then, system fl5.20p - 
f l5.22p is hyperbolic in the region where the angle 9 between VL and the direction of 
propagation is less than this limit value, independently of the density p. Fig. 15.21 summa- 
rizes the different the types of macroscopic limits of the system in dimension 2, when the 
density p, and the angle 9 between VL and the propagation direction vary. The behavior 
of the system at the crossings either between the hyperbolic and non-hyperbolic regions 
or between the ordered and disordered regions, remains an open problem. We note that 
non-hyperbolicity problems appear in other areas such as the motion a an elastic string 
on a plane [2H]- 



6 Conclusion 

In this paper, we have derived a macroscopic model for particles undergoing self-alignment 
interactions with phase transitions. This model is derived from a time-continuous version 
of the Vicsek model. We have identified two regimes. In the disordered regime, the 
macroscopic model is given by a nonlinear diffusion equation depending on the small 
parameter e describing the ratios of the microscopic to macroscopic length scales. In the 
ordered regime, the model is given by a hydrodynamic model for self-alignment interaction 
which is not hyperbolic. Many problems remain open. Among others, a first one is to 
determine the evolution of the boundary between the ordered and disordered regions and 
to understand how the models in the two regions are connected across this boundary. The 



20 



TT 




' ^ ' ' ' 

2 4 6 8 

Density p 

Figure 4: Types of macroscopic limits in dimension 2. Around the threshold value p = 2, 
none of the diffusion or hydrodynamic limit is valid. The study of this transition is still 
open. 



second one is to understand how to cope with the non-hyperbolicity of the model in the 
ordered region and possibly modify it by adding small diffusive corrections. Numerical 
simulations of the particle model are in progress to understand the behavior of the model 
in the two regimes. 

Appendix 1. Poincare constant 

In this appendix, we prove the following: 

Proposition 6.1. We have the following Poincare inequality, for if) G H^{B>): 

(|V.^Hm.o > A.((V^ - WA/«a)')A/..- (6.1) 

The best constant in this inequality is the smallest positive eigenvalue of the operator 

^:n = -T^V.-(M.nV.-)- (6.2) 

We define the linear operator L* hy 

Ll{g){e) = -(sin^)2-"e-"^°^^((sin^)"-V'=°^VW)'- (6-3) 
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Then one of the following three possibilities is true: 

(i) is the smallest eigenvalue of the Sturm- Liouville problem 

Llig) = Xg, (6.4) 

for g G C^([0,7r]) with Neumann boundary conditions (g'{0) = 5''(vr) = 0) and such 
that {sm6)'"~'^e'^^°^^ g{6)d6 = 0, and the eigenspace of L*^^ associated to the eigen- 
value is of dimension 1, spanned by u ^ h^ni^ ' where the function 9 ^ hQ[cos6) 
is smooth, positive for ^ 9 < 6q and negative for 6q < 6 ^ n . 

(a) Ak is the smallest eigenvalue of the Sturm- Liouville problem 

LUg) = LUg) + ^^g{e) = Xg, (6.5) 

forg G C^([0, vr]) with Dirichlet boundary conditions (g{0) = g^ir) = 0), and the eigenspace 
ofL*^Q associated to A^ is of dimension n—1, consisting in the functions of the form ip = 
h].{u: ■ VL)A ■ u for any vector A G M" such that Q ■ A = 0, with 9 i— t- h^{cos6) a smooth 
positive function for < 6 < n . 

(Hi) The two above Sturm- Liouville problems have the same smallest eigenvalue A^, and 
the eigenspace of L*^ associated to A^ is of dimension n, spanned by the two types of 
function of the above cases. 

Proof. First of all, we have 

(|V^^Hm.o > (minM^n) / |V^^p ^ (min M«f,) (n - 1) lU'lA , (6.6) 

and 

((^ - ^ ((^- / ^)')a/.„ ^ (maxM^f,) f L - U\ . (6.7) 

Js Js \ Js J 

The second inequality of fl6.6p follows from the Poincare inequality on the sphere: 

1 



Js J n-1 
The first inequality of (16. 7p follows from the fact that 

((^ - f ^f)M^a - ((^ - WKU,f)M^a ={1^-1 ^M,^ > 0. 



Eqs. (16. 6p and (16. 7p lead to the Poincare inequality (16. ip with 

min M«,o 

A, ^ n - 1 -— = n - 1 e'". 

max M«;Q 
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We use the inner product (</?, V^) {(pip)M^n^ adapted to M^^- We denote by L'liS) 
(resp. Hl{S)) the functions V e L'^{S) (resp. in such that (V^)A/«a = 0- 

The operator L*^ given by (16. 2p is self-adjoint since (V^j'?/' ■ Vaj<y5)j\/^n = m^q- 
It is then easy to see, using Lax-Milgram theorem, that if ip belongs to Ll^(S>) then there 
is a unique solution ip G Hl(E>) to the equation L*^q4' = V^- The so-obtained inverse 
operator is then compact and self-adjoint. By the spectral theorem, we get a basis of 
eigenfunctions, in the Hilbert space ^^^(S), which are also eigenfunctions of L*^^. If we 
denote the largest eigenvalue of the inverse of L*^, then it is easy to see that is 
the best constant for the following Poincare inequality, in the space Hj.{E>): 

Since the constants trivially satisfy this inequality, this shows that is the best constant 
for the Poincare inequality (16. ip in if^(§). 

The goal is now to reduce the computation of the eigenvalues to simpler problems, 
using separation of variables. We write cu = cos^f2 + sm9v, where v belongs to the unit 
sphere, orthogonal to VL. We identify Vt with the last element of an orthogonal basis of M*^, 
and we write v G Sn-2- 

By spherical harmonic decomposition in an adapted basis (see for example [TH], ap- 
pendix A), we have a unique decomposition of the form 

^(c^) = Y^atrntivl (6.8) 

k,m 

where {Z^{v))keii,km\ is a given orthonormal basis of the spherical harmonics of degree m 
on for m G N, with km = ("^^2^^) ^ ^^1^2'^)- ^ continuous, g'^ is given by 

g'^{e)= f ^{cos9n + smev)Z^{v)dv. (6.9) 

JSn-2 

We now show that the decomposition (16. 8p remains stable under the action of the opera- 
tor Lf^Q, so that its spectral decomposition can be performed independently for each term 
of the decomposition. 

First, we examine the case of dimension n ^ 3. Let tp{u}) = g{6)Z{v). We have 

V^^ico) = g\e)eeZ{v) + ^V,Z{v), 

sm 6 

where the unit vector eg is given by 

ee = = — - u ^ uj)n. 
sm (J 

We take functions ip^u) = gi{9)Z^{v) and ^{u) = J2k,nJMZt{v). Since the 
spherical harmonics are orthonormal, and are eigenfunctions of A^, for the eigenval- 
ues —m{rn + — 3), we get: 

(v.^ ■ v.^)m.. = HflL'mt'ie) + '^'^^^fje)gt{e)]{smer-\^-'de. (6.10) 
^0 
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Suppose m ^ 1. Then, it is easy to see that the function belongs to Hl{E>) if and 
only if {sm9)^^^g' G L^(0,7r) and (sin6')t~2^' ^ L^(0,7r). This condition is equivalent 
to the fact that g ^ V, where V is defined by (15. 5p . and which we denote by for 
convenience: 

= {g I {sm9)^^'g e L\0,n), (sin^)?"!^ e H',{0,n)}. 

Suppose now that m = 0. Then is a constant, and the condition ip G -f^^(S) 
is equivalent to the first condition only: {sin 6)^~^g' G L^(0,7r), under the constraint 
that J^{sme)''-^e''''°^^g{e)de = 0. We will denote this space by V^: 

= {^1 (sin^)t-V e L^{0,7r), J^{sme)''~^e^'°''^g{e)de = 0}. 
Formula fl6.10p then suggests to define the operator ^ : V^' — j- (V^)* by 

r f{9)Ll^,g{9){sm9r-'e^^"^'d9 = f [f g' + "'^lyV ^] (sin g)-^ e^^^ dg. 

(6.11) 

From (l6.1Up . it follows that, if we decompose ip{uj) = J2km9rni^)^mi'^)y then 



showing that L*^ is block diagonal on each of these spaces (tensorized by the spherical 
harmonics of degree m on Sn-2)- So we can perform the spectral decomposition of L*^ 
by means of the spectral decomposition of each of the L*^rn- is indeed easy to prove, 
using Lax-Milgram theorem, that the operators L* ^ have self-adjoint compact inverses 
for the dot product {f,g) = fg {sin 9)"'~'^e'^'^°^^d9 . Therefore the eigenf unctions and 
eigenvalues of L*^ correspond to those of the operators L*^rny m G N. If we denote 

by XK,m the smallest eigenvalue of L* ^, we finally get 

= min{A«;,m,"^ G N}. 

We notice that 



A«,™ = inf I £ f{9) LIJ{9) (sin 9)^"^ e^^^^' d9 



m 



f{9) (sin e^'°'^ d^ = 1 1 



but since all the V^^ are the same for m ^ 1, and since 

/ -7e/^(sin^)"~^e'^™'^Me ^ / (sin e)""^ e'^™'M^, 
Jo Jo 

we get 

A«:,m+i > K,m + {m + l){m + u - 2) - m{m + n - 3) = A^.m + 2m + n - 2. 
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Finally, is the minimum between o and A^,!- The eigenfunctions for the operator L*^q 
being smooth, this is also true for the operators L*^^, by formula fl6.9l) . So we can 
transform the definitions (16. lip by integration by parts. 

Indeed, if go is an eigenfunction (in V^) associated to L* q and an eigenvalue A, then go 
is smooth and satisfies the Sturm-Liouville eigenvalue problem 

Kgoid) = -(sin^)2-"e-''^°^^((sin^)'^-2e'=™^V(^))' = A^?o(^). 

Conversely, a smooth function with the condition JJ^(sin 6')"~^e^'^°'^^5'(^)d^ = belongs 
to V^. Actually, in dimension n ^ 3, we do not need to impose the Neumann boundary 
conditions: they appear naturally, since we have 

J^k9o - -e (e g^) - ^5'o - ^fl'o- 

Therefore by continuity at ^ = and vr, g'o^O) = g'oij) = 0. Then, using classical Sturm- 
Liouville oscillation theory (see |37] for example), we find that the first eigenspace of L* 
is of dimension 1, spanned by a function gKfl{d)i which is positive for ^ ^ < 6*0 and 
negative for 6o < ^ n. 

Similarly, if gi is an eigenfunction (in V^^) associated to L* ^ and an eigenvalue A, then gi 
is smooth, with gi{0) = 5'i(7r) = and satisfies the Sturm-Liouville eigenvalue problem 

LiMO) = L*Me) + ^,gm = ^9iio). 

And conversely, if a function satisfies Dirichlet boundary conditions while being in C^([0, tt]), 
it belongs to V^. Once again, if ^ 3, we do not need to impose the Dirichlet boundary 
conditions in the C^([0,7r]) framework, since we have 

Llg, = -e—\e-^%[)' - ^,g[ + ^^g, = Xg,. 

So, by continuity at 6* = and tt, 5'i(0) = (71 (vr) = 0, and then a first order expansion 
shows that continuity holds whatever the values of gQ{6) at the endpoints are. Again, using 
classical Sturm-Liouville theory, we find that the first eigenspace of L* is of dimension 1, 
spanned by a function 5'k,i(6'), which keeps the same sign on (0, vr). 

The case Ak,o < ^^,1 corresponds to case (i) of the proposition. Since a spherical 
harmonic of degree on the sphere Sn-2 is a constant, introducing /i^ such that (cos 6) = 
QKfiiG) allows us to state that the eigenspace of L*^^ associated to the lowest eigenvalue is 
spanned hy u ^ h^^iuj ■ Vt) . 

The case A^ > ^k,i corresponds to case (ii) of the proposition. The spherical harmon- 
ics of degree 1 on the sphere §„_2 are the functions of the form v ^ A - with A ■ f2 = 0. 
Introducing /i^ such that h^{cos6) sin 6 = gnfl{9) allows us to state that the eigenspace 
of L*Q associated to the lowest eigenvalue is of dimension n — 1, consisting of the functions 
of the form u ^ h],{(jj ■ VL) A ■ with A any vector in M" such that A ■ = 0. 

Finally, the case Ak,o = 1 corresponds to case (iii) of the proposition and this ends 
the proof in the case of dimension ^ 3. 
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We now examine the special case of dimension n = 2. We identify H^{E>) with the 27r- 
periodic functions in Hl^^(M.). So, is the smallest eigenvalue of the periodic Sturm- 
Liouville problem 

LUg) = L*M = _e— os.^g.cos = Xg, 

for functions g such that J^^e'^'^°^^g{6)d6 = 0. Here the decomposition corresponding 
to fl6.8p is the even-odd decomposition (there are only two spherical harmonics on Sq^ 
the constant function of degree and the odd function of degree 1). The odd part go 
of g can be identified with a function of Hq{0, tt), and it is easy to see that the odd part 
of L*^{g) is L*^{go), and similarly for the even part ge- So, we can perform the spectral 
decomposition of L* separately on the spaces of even and odd functions. 

Actually, if is a solution of the Sturm-Liouville periodic problem, the function g{6) = 
^-Kcoseg^g^^ — 6) is another solution with the same eigenvalue. Furthermore, if g is 
odd, then g is even and conversely. So the eigenvalues are the same for the odd and 
even spaces problems. Therefore, in dimension n = 2, proposition 16.11 can be refined 
and we can state that case (iii) is the only possibility: the eigenspace of L*^ associated 
to Ak is of dimension 2, spanned by an odd function g°, positive on (0,7r), and an even 
function g"^ = positive for < 6 < 6q and negative for 6q < 6 < tt. The proof of 
Proposition 16. II is complete. ■ 

We can now state a conjecture, which refines proposition 16.11 if true, and which is 
based on numerical experiments. 

Conjecture 6.1. (i) When n > and n ^ 3, only statement (ii) of Proposition \U7T\ is 
true. 

(ii) The function k ^ K,^ is increasing. 

We also observe numerically that Ai ~ k when kappa is large. 

Some investigations are in progress to prove the monotonicity of the eigenvalue with 
respect to k, based on formal expansions similar to those used in Section 5 of |18] . 



Remark 6.1. At the end of the proof of Proposition \ 6.1\ we have seen that in dimen- 
sion n = 2 only statement (iii) is true. The proof uses a transformation of the solution of 
an eigenvalue problem into the solution of another eigenvalue problem. We can try to find 
a similar transformation in dimensions n ^ 3: if f satisfies L^ of = Xf (with Neumann 
boundary conditions) then f = e~'^^°^^ dgf{n — 6) (with Dirichlet boundary conditions) 
satisfies 



/Li/(sine)"-2e'^™'^de = A / (sin^)"-^ e''™^'' d^ 

/»7r 

-K{n-2) / cose/2(sin^)"-2e'^^°'^^d^, 
Jo 

so if we can prove that cos9 f^ (sin 6*)""^ e'^'^°'^^ d6' > 0, we can deduce that Aq > Ai. So 
far we have been unable to prove this estimate. 
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Appendix 2. Numerical computations of the coeffi- 
cients 



We adopt a finite difference approach to compute the function g,^ associated to the GCI's 
and defined by 05.61) . We consider the function such that /k(^) = {sm9)2~^g^(^9). 
In particular, since & V defined by (15. 5p . belongs to Hq^O,^). Since g^, satisfies 
(15.61) . /k satisfies 



— e 



-K cos 9 f^K COS 9 fl \l 



;2^(i+¥^os^ 



K COS 6) fi^ = sin 2 6. 



We discretize the interval (0,7r) with + 1 points 9i = jjin, and denote by /* an 
approximation of at these points. Since G iJg(0,7r), = = 0. We define ej. = 
gKcosS,^ A second order approximation of (e'^™'^^/^)' at 9i is then given by 



7' 

J K 



Introducing 



dl 



n 



2 sin^ Oi 
hi = 



[1 + ^COS^^i) - KCOsOi + 



Ar2 



and bl 



Ar2 



the vector F = (/*)jg|[i^7v_i]] is the solution of the linear system AF 
vector S is (sins' 6'j)jgji jv„i|, and the tridiagonal matrix A is defined by 



S, where the 



A 















\ 


K 


dl 


bl 











bl 


dl 































2 











d^- 
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(6.1) 



We use the trapezoidal method to perform the integrations in the definitions (13.1 ip and 
(15. 9 p of c and c. The other coefficients p, A and 6^ are then directly computed from c 
and c. The numerical results provided in Figures IH1 I5.2I have been obtained for = 3000. 

We now detail how we obtain an approximation of the Poincare constant A^. By 
Appendix 1, A^ is the minimum between Ak,i and Ak,o, which are the smallest eigenvalue 
of two Sturm-Liouville problems. Several algorithms exist to compute eigenvalues of 
singular Sturm-Liouville problems (which is the case here whenever n ^ 3) with a good 
precision [3]. However, we use a simpler method based on finite differences. 
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Actually, A^,! is the smallest eigenvalue associated to problem fl6.5p . with g E V. So, 
considering once again the function / such that f{9) = {sin9)^~^g{9), the vector AF, 
with A defined by f l6.ip . gives a second order approximation of {sin 9)^~^L*^g (9) = Xf{9) 
at the points 9i. So we can take the smallest eigenvalue of A as an approximation of A^^i. 

We now look for an approximation of X^^. Let g he a. solution of the Sturm-Liouville 
problem f l6.4p with Neumann boundary conditions. We introduce G = )jg|[o,7v-i]. 



the vector of approximations of g at the points 9^,i 



N ' 



vr. 



Introducing mj. 



(sin 6*)" ^ gKcos6»,^ ^ second order approximation of L* c/ at the point 6'j_,_i , with i G |1, iV— 2] 
is then given by 



i+i 



^+1 



With the Neumann boundary conditions, the approximations at the points 9i and 9j^_ 
are given by 



Introducing 



AT- 



TV- 



^2 



ml 



i+4 



and 



" 2 



a second order approximation of L*^g is given by BG, where the tridiagonal matrix B is 
defined by 



B 



l-bl 


1 







3 


3 


hi 




&I 




















V 







Ok 



N- 
~N- 



(6.2) 
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So we can take the smallest positive eigenvalue of B as an approximation of X^fl 
(excluding the constant functions). The computations of Fig. [T] have been performed 
with = 300 points. 
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